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Abstract 



The data on average hadron multiplicities in central A+A collisions measured 
at CERN SPS are analysed with the ideal hadron gas model. It is shown that the 
full chemical equilibrium version of the model fails to describe the experimental 
results. The agreement of the data with the off-equilibrium version allowing for 
partial strangeness saturation is significantly better. The freeze-out temperature of 
about 180 MeV seems to be independent of the system size (from S+S to Pb+Pb) 
and in agreement with that extracted in e + + e~, p + p and p + p collisions. The 
strangeness suppression is discussed at both hadron and valence quark level. It is 
found that the hadronic strangeness saturation factor 73 increases from about 0.45 
for p + p interactions to about 0.7 for central A+A collisions with no significant 
change from S+S to Pb+Pb collisions. The quark strangeness suppression factor 
Ag is found to be about 0.2 for elementary collisions and about 0.4 for heavy ion 
collisions independently of collision energy and type of colliding system. 
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1 Introduction 



Many years of experimental effort in the field of high energy nuclear collisions yielded a 
large amount of data on particle production at different collision energies (up to 200 A 
GeV/c) and for different colliding systems [[J. These results allow studying the proper- 
ties of strongly interacting matter at high energy densities. Ultimately, at high enough 
collision energy one expects to create in the laboratory the Quark-Gluon Plasma (QGP), 
a form of matter in which effective degrees of freedom are quarks and gluons instead of 
hadrons and hadronic resonances. The formation of QGP with deconfinement of quarks 
and gluons should hopefully be reflected in the final hadron production, provided that 
the expected modifications of the entropy and strangeness content of the system survive 
hadronization and reinteractions between final state hadrons. In principle the system evo- 
lution is determined by QCD. Nevertheless, the formation of hadrons is a process entirely 
lying in its non-perturbative domain, hence, in order to study the final state, one has to 
resort to phenomenological models such as string or statistical (thermal) models. 

The statistical models, whose prototypes 0, [|, || date back to '50s and '60s, are based 
on the assumption of local filling of available phase space according to statistical laws, once 
collective effects have been taken into account. This ansatz allows the characterization of 
hadron production by means of few parameters such as temperature, volume and chemical 
potentials. Furthermore, parameters accounting for possible departures from complete 
equilibrium are often introduced. 

A long and rich history of thermal models is related to their surprising success in the 
description of many aspects of high energy collisions || . It has been shown recently that 
the ideal hadron gas model allowing for non-equilibrium strangeness abundance is able 
to reproduce hadron multiplicities in e + + e~ , p + p and p + p interactions over a large 
collision energy range |§, 0, §. In this paper we apply the same model to analyse the 
data of hadron multiplicities in central nucleus-nucleus collisions at CERN SPS collision 
energies. The data on Sulphur-nucleus collisions at 200 A GeV/c were already analysed 
using thermal models in many previous works (for a review see Ref. 0). We add to 
this data set the preliminary results on hadron production in central Pb+Pb collisions 
at 158 A GeV/c. The use of the same formulation of the thermal model for the analysis 
of data from p + p to central Pb+Pb collisions allows us to study the evolution of the 
parameters of the model in the full range of the colliding systems. Our analysis is based 
on hadron multiplicities integrated over the full phase space because their use minimizes 
the influence of collective effects on final results allowing a formulation of the model to 
be tested with a minimal number of assumptions and parameters. 

The paper is organized as follows: in Sect. 2 we describe the hadron gas model which 
is confronted with the experimental data in Sect. 3. Discussion and conclusions are given 
in Sect. 4. 

2 Hadron Gas Model 

In this section the hadron gas model used in the present analysis is sketched; a more 
detailed description can be found in Refs. |7], |[. The model postulates the formation of 
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an arbitrary number of hadron gas fireballs each having a definite collective momentum as 
the result of an interaction between the two colliding systems. The parameters describing a 
hadron gas fireball at thermal and chemical equilibrium are the temperature 7$, the volume 
Vi in its rest frame, and the quantum numbers, i.e. electric charge Qi, baryon number 
Bi and strangeness Si. Charm and beauty have been excluded from the calculations as 
the thermal production of heavy-flavoured hadrons is negligible with respect to non- 
heavy-flavoured hadrons in the expected range of temperatures, namely 100-i-200 MeV. 
If Q° = (Qi, Bi, Si) is the vector of i th fireball's quantum numbers and N is the number 
of fireballs, the following constraint: 

N 

EQ° = Q°> (i) 

i=i 

where Q° is the quantum vector fixed by the initial state, must be fulfilled. The average 
yield of any hadron in the i th fireball can be derived from the partition function: 

Zi(<$) = E e- Ei/Ti 6 Qi ,Q9 , (2) 

states 

which is calculated in the canonical approach, i.e. by using only the multi-hadronic states 
having the same quantum numbers of the fireball. 

A non-equilibrium parameter 734 accounting for a possibly incomplete strangeness 
chemical equilibration is introduced by multiplying by 7^ the Boltzmann factors e~ Ej ' Ti 
associated to the j th hadron in the partition function, where s is the number of its valence 
strange quarks and antiquarks. Although this factor was introduced heuristically [[UJ 
and used as a purely phenomenological parameter in the analysis of elementary collisions 
||, [7], H, it can be shown that 73 formally is the fugacity associated to the number of 
strange + antistrange quarks in the hadron phase in a grand-canonical framework . 

The overall average multiplicity of each hadron species is the sum of all average yields 
in each fireball. In principle any configuration {Q?, . . . , Q%} of fireballs in the event may 
occur, so that average hadron abundances depend on the probability w(Q^, . . . , Q° N ) of 
occurrence of a given configuration besides the whole set of fireball thermal parameters 
Ti, Vi and 7sj. However, it can be shown that if such configuration weights w(Q?, • • • , Q° N ) 
are chosen in a statistical fashion |7|, ||: 

w(O 0° ) - ^Q°.Q° II=i^(Q°) 

LQ0,...,Q$£ a Q°,Q0 ll i= l Ai[\4i ) 

and the fireball freeze-out temperatures and 7si suppression factors are the same, namely 
Ti = T and 7s« = 7s, then the average hadron abundances nj at freeze-out depend only 
on the global volume V = Vi and Q° through the following equation: 

where the upper sign is for fermions and the lower for bosons; the function Z is the global 
partition function f7|, § and is the quantum number vector of the j th hadron species. 
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The special choice of weights w(Q?, . . . , Q%) in Refs. 0, || leads to the same expres- 
sion of average multiplicities relevant to a system in global equilibrium even if the fireballs 
are not in mechanical equilibrium. It might be argued that the difference between the 
rapidity spectra of baryons and antibaryons - existing even in central A+A collisions J[^] 
- question the validity of this choice. Nevertheless this choice has a remarkable property, 
namely it removes the dependence of hadron average multiplicities on both the number 
of fireballs and their ordering either in size or in space (reabsorbed in the global volume) . 
As a consequence, much freedom is left to possibly reproduce the rapidity distributions 
within the model keeping the same quantitative expressions for hadron abundances (see 
Appendix A). However, even if the choice of the weights w(Q?, . . . , Q%) was not cor- 
rect, the corrections to eq. @ are expected to be small because the relative abundances 
are predominantly determined by the intensive thermal parameters and not by fireball 
quantum configuration weighting. 

The chemical factors Z(Q,° — Zqj)/Z(Q°) in eq. (^) implement the dependence of the 
yields upon the chemistry of the system and replace the chemical potentials in the proper 
canonical formalism. For very large volumes, as expected in heavy ion collisions, it can 
be proved || by means of a saddle-point approximation, that eq. (|) reduces to: 

n, = (2Jj + 1) g£(Tl) m 7s^K 2 (^) e^A- 1 */* e -^A-^/4 j (5) 



"j 



where A is a 3 x 3 matrix whose elements are proportional to V: 



1 ^(24 + 1) I 3 j' s 'e-V^' T 

where the sum runs over all hadron species. In eq. @ the chemical factors are transformed 
into a product of two factors: the first one can be written as exp[//x • cy/T] where \x is 
a traditional set of chemical potentials, whereas the factor exp[— / 2 q ? A _1 q J /4] has no 
grand-canonical corresponding quantity. Its presence entails a suppression of hadrons 
having non-vanishing quantum numbers with respect to neutral ones, owing to the finite 
size of the system. Indeed this factor takes its origin from the requirement of exact 
conservation of initial quantum numbers. In the infinite volume limit A -1 goes to zero 
and the grand-canonical formalism is fully recovered. 

An important problem to face in modelling heavy ion collisions is the fact that particle 
multiplicities are measured by averaging over events with a varying number of participant 
nucleons. For central collisions of identical nuclei (S+S and Pb+Pb in this paper) the 
fluctuations in the number of spectator nucleons, Aspec, can be considered to follow a 



Poisson distribution fll3[ , well approximated by a Gaussian distribution if (Aspec) 3> 1- 
Hence fluctuations in the number of participant nucleons, Apart = 2 ■ A — Aspec, can be 
described by a Gaussian distribution with mean value (Apart) and variance (A SPEC ). For 
the considered collisions (S+S and Pb+PbQ) it is (A SPEC ) <C (Apart), implying that the 
mean of participant nucleons distribution is large in comparison with its width. Assuming 

4 Validity of the above consideration for central S+Ag collisions is questionable as fluctuations are 
probably dominated by a geometrical effect. 
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that T and 7s do not depend on the number of participant nucleons, the average number 
of directly produced hadrons of species j is: 



(n,) = (2J j + 1) W dVF(Q°,V) ^E(Tl) m 7? J ^K 2 (^) e^*" 1 ^ 4 , 

QO J Z7F i=l ' 1 



(7) 

where -F(Q°, V) is the joint probability of observing an event with global volume V 
and quantum vector Q° (determined by the number of protons and neutrons involved in 
the collision). The removal of the integral from eq. (6) by introducing mean chemical 
potentials JT and mean volume V must be undertaken with great care since such mean 
values in general would depend on the hadron species and would not be the same for all of 
them. This is clearly understood because /xand A, depending on Q° and V, multiply q,, 
in a non-factorisable form. However, if a reasonable ansatz of a narrow Gaussian shape 
for the distribution function F(Q°, V) is assumed, according to previous discussion, it 
can be shown (see Appendix B) that the simple averaging procedure introducing mean 
chemical potentials and volume can be used provided that A -1 <C 1. In other words, for 
the following simple averaging formula to hold: 

(n,) « (2Jj + 1) YL 7s^K 2 fe e^/* , (g) 

a nearly grand-canonical regime and reasonably small fluctuations of V and Q° are needed. 
How small such fluctuations must be, will be checked a posteriori in the actually examined 
collisions. The hadron average multiplicities are calculated with eq. (§) in which V, T, 
7s and ~p, B are free parameters to be fitted to the data, while ~fi s and JIq, the strangeness 
and the electrical chemical potentials respectively, have been determined by means of 
two additional constraints: the strangeness neutrality and the conservation of the ratio 
[Z\ + Z2) j {A\ + A2) formed with the atomic and mass numbers of the two colliding nuclei: 

j 

It should be noted that the first of the two constraints in eq. @ is valid on an event by 
event basis whereas the second, only averaged over a large number of collisions. 

All previous equations are concerned with primary hadrons, namely particles and 
resonances directly emitted from the hadronic source and not coming from secondary 
decays. On the other hand, since actual measurements include feeding from heavier 
hadrons and resonances, the hadron production rates to be compared with the data have 
been calculated by letting all primary hadrons decay according to known branching ratios 
until particles considered stable by the experiments are reached. Among primary hadrons 
we included all particles and resonances up to a mass of 1.7 GeV; the masses of resonances 
have been distributed according to a relativistic Breit-Wigner. The needed values of 
hadron masses, widths and branching ratios have been taken from the most recent Particle 
Data Book WA . 
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3 Comparison with the Data 



3.1 Full Equilibrium Model 

The hadron gas model described in the previous section has been used to fit the data 
on hadron abundances in central S+S, S+Ag at 200 A GeV/c and central Pb+Pb colli- 
sions at 158 A GeV/c measured by NA35 and NA49 Collaborations at CERN SPS. We 
used average hadron multiplicities measured in (or extrapolated to) full phase space. The 
compiled data and the references to the original papers, where the experimental details 
(acceptances, extrapolation procedures) can be found, are given in Table |l|. Since pre- 
liminary data on central Pb+Pb collisions are still poor, we decided in this case to use in 
addition two particle ratios measured in the central rapidity region, namely K + /K~ and 
A/A. The measured rapidity distributions of these particles in the acceptance region are 
similar to those of the corresponding antiparticles []2T|, thus justifying our decision to use 



those ratios as estimates of the full phase space ones. 

The analysis of particle multiplicities in e + + e~, p + p and p + p collisions within a 
hadron gas model indicates that the strangeness production in such elementary collisions 
is not high enough to ensure the complete local chemical equilibrium at hadron level. It 
may be expected that in the case of central A+A collisions, due to a much larger volume 
of the interaction region and an increased role of hadron rescattering, a full local chemical 
equilibration at hadronic level can be attained. 

In order to test this hypothesis we started the comparison of the experimental data 
with the hadron gas model by using first its fully equilibrated version, i.e. by setting 
7s = 1. We performed a graphical test similar to that made in Ref. [22| by plotting in the 



T-~p B /T plane the bands determined by the central values of some of the most relevant 
ratios of hadron yields and by their lcr variations. As the overall multiplicities are very 
large, we set A -1 = in eq. (^); in this grand-canonical limit the multiplicity ratios 
depend only on the intensive free parameters T and ~p. B /T as the mean volume V cancels 
out. This approximation turns out to be satisfactory for all the examined collisions as it 
is demonstrated below. In Fig. [I] a-c such bands are shown for S+S, S+Ag and Pb+Pb 
data, respectively. There is no evident common crossing region for all bands, indicating 
absence of complete equilibrium even in central collisions of nuclei as heavy as Lead. It 
should be noted that a crossing region does exist for Pb+Pb collisions with T ~ 120 
MeV and ~p B /T ~ 2.7, as long as only full phase space multiplicities are considered, i.e. 
excluding the A/A band. However, these latter values imply an antibaryon/baryon ratio 
of less than 10 -2 at primary level, which is unrealistically small with respect to the same 
ratio in S+S and S+Ag collisions. In fact the measured A/A band is quite far from the 
crossing point. 

In order to confirm the previous finding we performed a least-squared fit to the data 
with T, V and ~Ji B jT as free parameters, using the canonical corrections and keeping 
7s = 1 fixed. The details of the fitting procedure are described in Appendix C. The 
results are shown in the first part of Table |2] and the fourth column of Table |l|. The 
X 2 /NDF's are about 6-j-8 with large discrepancies between fitted and measured values. 
Therefore we conclude that a full equilibrium version of the ideal hadron gas model fails 
to reproduce full phase space hadron multiplicities in central A+A collisions at CERN 
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3|. This discrepancy can 



SPS energies. 

Our conclusion is in contradiction with the findings in Ref. 
be explained by the fact that the analysis in Ref. [23] was performed by using particle ratios 
in various regions of rapidity and transverse momentum which, unlike in our analysis, 
requires additional dynamical input beyond a simple statistical ansatz. Secondly, a proper 
statistical comparison between model predictions and data has not been performed in 
Ref. ]|23 



3.2 Off-Equilibrium Model 

We tested the off-equilibrium version of the ideal hadron gas model repeating the fits 
with 7s as a free parameter. The results of these new fits are shown in Table [1] and in 
the lower half of Table |2| while the comparison between fitted and measured multiplicities 
(or ratios) is shown in Table |1| There is a good agreement for all particles with some 
exceptions such as antiprotons in both S+S and S+Ag collisions. The x 2 /NDF for S+S 
and S+Ag collisions is about 3.5, thus significantly lower than for the full-equilibrium fit. 
The x 2 /NDF for central Pb+Pb collisions is close to 1. 

The obtained temperatures and baryon chemical potentials are quite compatible with 
a common value and so are the 7s values which turn out to be definitely less than 1 in all 
the three examined collisions. Also quoted in Table ^| are the obtained chemical potentials 
~p s and ~pq and the range of variation of the matrix A -1 elements; their smallness bears 
out the saddle-point approximation used in this analysis and the previously described 
graphical test by confirming the proximity to the grand-canonical regime. 



3.3 Discussion of the Analysis 

The dependence of the fitted parameters on the hadron mass spectrum cut-off has been 
checked by repeating the fit with lower cut-off values and found to be negligible. 

The validity of the approximated averaged formula in eq. (§) in the presence of par- 
ticipant fluctuations has been checked by calculating its first order corrections according 
to the formulae quoted in Appendix B with the assumption of moderate Gaussian fluctu- 
ations of baryon number B, electric charge Q and volume V. The corrections have been 
estimated repeatedly by random variations of the correlations, assumed to be positive, 
between B, Q and V in order to find out a maximum value. Particles undergoing the 
most significant variations with respect to the average yield are baryons whose production 
increases mainly owing to terms proportional to (5V/V) where 5V is the volume disper- 
sion. Such variations turn out to be almost constant for the three kinds of collisions as a 
function of the relative dispersions of B, Q and V. They are within 5% with no volume 
fluctuation even for 5B/B and 5Q/Q of 30% but can raise to more than 30% if also a 
volume dispersion SV/V = 30% is included. 

As follows from the previous discussion (see Sect. 2) and the data presented in Table [1|, 
the relative width of baryon number distribution is expected to be around 6% and 2% 
for central S+S and Pb+Pb collisions, respectively. For these participant fluctuations 
the corrections to baryon yields range from about 2.5% (1%) for SV/V = 10% up to 
17% (13%) for SV/V = 30% in S+S (Pb+Pb) collisions. It is worth remarking that the 
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B, Q and V fluctuations would mainly influence the determination of the baryochemical 
potential in the fit as antibaryons and other particles production are much less affected 
by them compared with baryons. 

Whilst S+S and Pb+Pb collisions involve two identical nuclei, S+Ag is an asym- 
metric nuclear collision and the use of the eq. (|9|) to determine strange and electrical 
chemical potentials may be not appropriate. According to eq. (§) the average ratio Q/B 
of participant nucleons in S+Ag collisions is about 0.45 whereas it is actually likely to 
be somewhat higher owing to the fact that Sulphur (Z/A—0.5) is the smaller nucleus. In 
order to prove that our results are independent of the previous assumption, we repeated 
the fit by varying the charge/baryon number ratio and taking the extreme case Q/B=0.5. 
We found T = 177.5 ±7.9 MeV, VT 3 exp(-0.7GeV/T) = 6.37 ±0.47, 7 S = 0.707 ±0.063, 
■p B /T = 1.325±0.081 with a x 2 /NDF = 7.0/2, a fit result not significantly different from 
that quoted in Table |j. 

A further technical problem in the x 2 fit is concerned with data redundancy. Whenever 
two or more data points are in a relationship that is not dependent on any free param- 
eter, then the x 2 significance might be questionable as the effective number of degrees 
of freedom is overestimated. This is the case for S±S collisions where the (approximate) 
relation (K°) = ((K + ) ± (K~))/2, owing to isospin symmetry and independent of model 
parameters, links three data points. However, it can be shown (see Appendix D), that 
such redundancy does not affect the extraction of thermal parameters nor their errors; its 
main effect is a lowering of the x 2 /NDF. 

Our analysis is done in the ideal pointlike hadron gas framework. Since the extracted 
temperatures are very high, corrections due to particle repulsion in principle should be 
considered. However, most corrections proposed in literature leave particle ratios un- 
changed so that the parameters T, 73 and chemical potentials are unaffected. On the 
other hand, our volume parameter V is the pointlike particle volume and should not be 
confused with the actual volume. In a recent analysis ]24| different particle ratios were 
compared with the full-equilibrium hadron gas predictions by introducing different hard- 
core radii for pions and other hadrons leading to an effective pion chemical potential. Such 
a procedure restores the agreement of the pion abundance which is always underestimated 
using thermal parameters extracted from the remaining hadrons |25j in a full-equilibrium 
model. However, the price to be payed is the introduction of additional, particle-species 
dependent, parameters. 



Our numerical results are quite different from a previous analysis of S±S collisions [p5 



which found T = 205 ± 39 MeV, 7 S = 0.95 ± 0.20 and fi B /T = 1.30 ± 0.17 by excluding 
negatively charged hadrons from the analysed data set. The fact that x 2 '/NDF > 3 in 
S±S collisions is expected to produce fluctuations larger than la in the fitted parameters 
if some data points are excluded in turn from the fit. In fact, by excluding negatives we 
found T = 200±10 MeV, 7s = 0.99±0.11 and fi B /T = 1.35±0.09, quite in agreement with 
Ref. |p5] . The remaining difference can be explained by the use of an updated data sample 



and particularly by the inclusion, in our analysis, of resonance widths; this especially 
enhances the p meson production and, consequently, the pion production improving the 
agreement with the measured negatively charged hadron multiplicity. 

Another recent analysis |26| of hadron abundances in Pb±Pb collisions found 7s = 
0.9 ± 0.09 but with the use of only strange baryons ratios in a limited phase space region 
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and omitting the multiplicity. 



4 Discussion and Conclusions 

The data on hadron multiplicities in central A+A collisions at SPS energies can not 
be reproduced by an ideal hadron gas assuming complete chemical equilibrium. This 
statement is also true for the central collisions of the heaviest nuclei, Pb+Pb. The hadron 
gas model supplemented with partial strangeness saturation agrees significantly better 
with the data as the ^/NDF reduces from about 7 to about 3 for S+S and S+Ag 
collisions and about 1 for Pb+Pb collisions. 

The behaviour of the intensive parameters obtained within the off-equilibrium model is 
shown as a function of the system size in Fig. [| where the result for the p + p data at 
similar collision energy || is also included. 

The lack of complete strangeness equilibrium at hadron level for central A+A collisions 
can not be interpreted as an effect of the choice of weights w(Q?, • • • , Qjv) described in 
Sect. 2 which is crucial to reduce the number of free parameters. Since 7s turns out to 
be < 1, it might be argued indeed that if the hadron gas fireballs were small enough and 
all with zero strangeness (so that the weights w(Ql, . . . , Q° N ) would be no longer those 
chosen in Sect. 2), a suitable canonical suppression || H could be generated without the 
need of 7s and a hadron gas in full chemical equilibrium would be recovered. Nevertheless 
this mechanism would have no effect on the yield of meson which is completely neutral, 
thus not suppressed by quantum number conservation at hadron level and having no 
known feeding from heavier light-flavoured resonances. Therefore, the measurement of 
(f) production in Pb+Pb collisions establishes the necessity of a significant strangeness 
suppression at hadronic level independently of the validity of the assumed fireball quantum 
configurations occurrence probabilities. 

An important conclusion can be drawn from the resulting chemical freeze-out temper- 
ature, which seems to be independent of the system size and it is, within errors, much the 
same as that in e + + e~, p + p and p + p collisions J7], ^|, ^7]]. This seems to indicate that 
the chemical freeze-out occurs close to the hadronization point and that the same mech- 
anism of statistical hadronic phase space filling at critical parameters of the prehadronic 
matter invoked as a natural explanation of elementary collisions results 0, [27|] also holds 
for heavy ion collisions. Moreover, the similar //b/T values for all studied A+A collisions 
suggest a common hadronization nuclear density. 

As the fit is not perfect in S+S and S+Ag collisions even in the off-equilibrium model, 
one may speculate that the small deviations from model predictions are due to secondary 
inelastic interactions between hadrons following the hadronization stage. While thought 
to be absent in p + p collisions, they may likely occur in A+A collisions where they 
can destroy the statistical character of the hadronization process. In fact, as inelastic 
cross sections for different processes are significantly different, hadron rescatterings may 
lead to decoupling of different particle species at different temperatures, thus affecting the 
single temperature fit. This mechanism is particularly well-suited to explain the observed 
deviation of antiprotons which may quickly annihilate in the baryon-dense medium formed 
in an A+A collision. On the other hand it should be mentioned that the observed small 
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deviations from the off-equilibrium version of the hadron gas model could simply stem 
from errors in the extrapolation procedures, from participant and volume fluctuations (as 
described in Sect. 3) and from the choice of weights w(Q?, . . . , Q%) (see Sect. 2). 

Whilst the temperature is constant, the strangeness suppression factor increases from 
about 0.45 for p + p interactions to about 0.7 for central S+S collisions at comparable 
nucleon-nucleon centre of mass energies. No further increase of the strangeness suppres- 
sion factor is observed for central Pb+Pb collisions. This observation has two important 
consequences: firstly, a heavy ion collision is not the result of an incoherent sum of nucleon 
collisions as far as strangeness production is concerned. In fact, due to isospin symmetry, 
7s must be the same in p + p and n + n collisions at the same yfs\ strangeness production 
in p + n interactions was measured to be the same as in p + p interactions ||28|| , thus the 



neutron content of colliding nuclei cannot account for the increase of 7s with respect to 
p + p and p + p collisions. Secondly, as canonical strangeness suppression was taken into 
account in extracting 73 in Ref. || and in the present analysis, the strangeness enhance- 
ment in heavy ion collisions cannot be fully attributed to the increased system size at 
hadron level. 

The relative production of strangeness has been intensively studied in elementary |2E 



and nuclear collisions [gO], £T| . It is usually expressed in terms of a strangeness suppression 
factor As defined as: 

(ss) 

S 0.5«uu) + (dd)) ' 1 ' 

where (ss), (uu) and (dd) are the mean multiplicities of newly produced valence quark- 
antiquark pairs at primary hadron level, before resonance decays. Thus the initial colliding 
valence quarks are excluded in calculating As- A major problem in the experimental 
determination of As is to account for unmeasured hadron abundances. The statistical 
model used in this analysis is a useful tool for this purpose because it reproduces well all 
measured hadron abundances both in elementary and nuclear collisions, thus providing 
a reliable quark counting method. In Fig. |3] As obtained by using model predictions for 
primary hadron multiplicities in e + + e~, p + p , P + P collisions with the parameters 
quoted in Ref. |27| and for A+A collisions is shown. It should be mentioned that in 
e + + e~ collisions the leading strange quarks in e + e~ — > ss events have been subtracted 
from the numerator of eq. (|H]) so that As contains only valence quarks created during 
the hadronization process. The As values for elementary collisions are consistent with 
a constant value of about 0.2, even for very high energy p + p collisions. The difference 
between 73 in e + + e~(~ 0.7) compared to p + p , p + p collisions (0.46-^0.56) as resulting 
from the off-equilibrium hadron gas model fit [[5], [8| is mainly due to two effects. As far 
as p + p collisions are concerned, the presence of six initial u, d quarks to be hadronized 
along with those newly produced brings about a 73 decrease. In fact, for constant T, V 
and 7s, a hadron gas with increasing baryon number and electric charge, i.e. increasing 
number of initial protons, has an increasing As (see Fig. |]). The physical reason is the 
lower energy threshold for strange pair production in a baryon rich environment where the 
dominant process is via N + it — > A + K while in the baryon free case it proceeds via kaon 
pair production. Secondly, T and 73 are anticorrelated in e + + e~ collisions. As the central 
fitted T value is lower in e + + e~ collisions in comparison with p + p and p + p collisions 
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2~Tfl , the central 7s value is expected to be higher in e + + e collisions to reproduce the 



measured strange hadron multiplicities. In fact, by repeating the same fit as in Ref. 
for e + + e~ collisions at y/s = 29 and 91.2 GeV, and keeping fixed T = 170 MeV instead of 
160.3 MeV and 163.4 MeV respectively, 7s turns out to be 0.69 and 0.62 instead of 0.72 
and 0.67 respectively with a slightly worse \ 2 - 

Our extracted value Ag is in agreement with a previous estimate based on quark 
counting method quoted in Ref. |3^] only for energies y/s < 100 GeV. On the other hand, 
the rise of As in p + p collisions claimed in Refs. [32, 33] is not observed. The reason of 
this discrepancy is the fact that As was estimated in Ref. |33| by using only K/tt ratio 
as experimental input and two parametrizations of hadron multiplicities [Q, |35[ which, 
unlike our parametrization J^, ||, do not satisfactorily reproduce all available measured 
multiplicities in p + p collisions Q. 

To summarize, the common characteristics of elementary interactions seems to be the 
independence of As on the collision energy, in the range examined in Refs. [|7j, |8| and 
type of colliding particles. This universal behaviour is broken in central A+A collisions. 
The value of As turns out to be about a factor two larger than the corresponding value 
for elementary interactions. Our value for A+A collisions is consistent with previous 



estimates based on quark counting method [p0| . Note that in Ref. |30| As was also 
estimated for p+A collisions and found to be consistent with that of p + p collisions 
and independent of the size of the target nucleus. This leads to the conclusion that no 
strangeness enhancement is observed in p+A collisions. 

The saturation of 7s and As factors as a function of the colliding system size for cen- 
tral A+A collisions suggests that the strangeness enhancement with respect to elementary 
collisions may already occur in the prehadronic phase and that secondary hadron scatter- 
ings, expected to be much more abundant in Pb+Pb collisions, are of minor importance 
for strangeness production. The strangeness enhancement effect observed in central A+A 
collisions and its independence of the colliding system size has been interpreted as due to 
Quark-Gluon Plasma formation in the early stage of the collision [[31], |36j already in S+S 
collisions and not only in central Pb+Pb collisions according to the interpretation of J j ip 
suppression J571 [58]. 
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Appendix 



A Rapidity Distributions and Statistical Weights 

In order to show that forward-backward peaked rapidity distributions for baryons and 
centrally peaked for antibaryons are not inconsistent with the assumption of statistical 
weights of eq. (|3|) we consider a simple example of p + p collisions. Since the derived 
expression of average multiplicities in eqs. does not depend either on the number 

of fireballs or their particular volumes Vi, . . . , Vn, we consider a toy model with three 
fireballs with equal volume Vf and sorted by boost velocities fa > fa > fa. According to 
the statistical choice of weights u?(Q?, . . . , Q%), owing to the equality of all parameters Vi, 
Ti and 7si of the fireballs, the probability of occurrence of baryon number configurations 
{1,0,1}, {0, 1,1} and {1, 1, 0} is equal; the same holds for more complex and less probable 
sets of configurations such as {2, 0, 0}, {0, 0, 2} and {0, 2, 0}. Therefore, as far as average 
hadron multiplicities are concerned, nothing changes if one replaces {0, 1, 1} and {1, 1, 0} 
with {1, 0, 1} and {0, 2, 0} with {2, 0, 0} in a half event sample and with {0, 0, 2} in the 
remaining half. The hadron abundances do not vary and a strongly forward-backward 
peaked rapidity distribution for baryons can be obtained as the fireballs having a non- 
vanishing baryon number, are always those in the forward or backward directions. 
In general this argument can be repeated for iV fireballs having an equal rest frame volume 
and an arbitrary set of ordered boost velocities fa > . . . > (3^. In this case the weights 
in eq. @ are symmetric: 

w(Ql {1) ,...,Q° a{N) ) = w(Q l ,...,Q° N ) (11) 

for any permutation a of the integers 1, . . . , N. Therefore, if p(Q?, . . . , Q° N ) are the actual 
weights, for eq. (H) to be valid, the condition to be fulfilled is: 

MQ?, • • • , Qat] = ^£KQ°(i), • • • , Q° ( ao) , (12) 

where the square brackets mean that the set [Q°, . . . , Q^] is a not-ordered one. This 
condition is weaker than a strict equality between w(Q5, . . . , Q^) and p(Q°, . . . , Qjv). 

To summarize, the compatibility between the expression for hadron multiplicities 
(eqs. (§,§)) and rapidity distributions can be achieved by choosing a model in which 
all fireballs have the same volume. Their boost velocities and their total number are 
allowed to vary event by event and can be determined by using actual hadron spectra. 

B Participant Nucleons and Volume Fluctuations 

In this section we point out the conditions to be fulfilled for the replacement of eq. (|7|) 
with its averaged version eq. (§) in the presence of fluctuations of participant nucleons. 
In general, the variation of the number of participants imply fluctuations of total baryon 
number B, electric charge Q and also global volume V of the colliding system. We 
assume that the associated distribution function F(Q, B, V) is a Gaussian and that the 
mean values Q and B are large, which is the case in the examined collisions. If the latter 
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condition is met, the sum over quantum vectors Q° in eq. (Utj can be turned into an 
integration: 



<"i> = (2^ + l)^£(Tl) m 7s S ^K 2 (^) 



x 



J dQdBdV F(Q, B, V) V e lfl ^ /T e^ A . (13) 



No more factor can be drawn out of the integral as the chemical potentials ^depend on the 
integration variables because of the quantum number conservation constraint J2j °Lj( n j) = 
Q° and the matrix A is proportional to the volume (see eq. @). Nevertheless, they can be 
expanded off the mean values Q, B and V up to first order, provided that the dispersions 
are not too large: 



~ M (X)+J M (X-X) 

Fih- 1 _ / V — V\ 

A- 1 ~ A-i + -^{V - V) = A-i (l - -^j . (14) 

where X = (Q, B, V) and J is the Jacobian matrix. Using the above expansions in eq. (|i~3|) 
one obtains: 



\ n 3 



00 r ( V — V\ 

) = E(4/dQdBdF + F(Q,B,V) 



x exp 



Iqj • J M (X - X)/T + 1% ^ V V V \ 3 /A 



(15) 



where {fij)i is just the I th term of the series in eq. The second term in the exponential 
is negligible if we are close to the grand-canonical regime and if the temperature is low 
enough to quickly suppress the terms of the series with high I: this condition is met for 
all hadrons if T < 200 MeV. Therefore: 



\ U 3 



00 r f XT — V\ 

)~J2(n j ) l JdQdBdV (^1 H = — J F(Q,B,V) expf/q, • J M (X - X)/T] . (16) 



If F is a multivariate Gaussian: 

1 



F(X) = ; exp[-(X - X) ■ C-^X - X)/2] , (17) 
(27r) 3 detC 



then the integral in eq. ( jig) can be solved analytically if the integration is extended to 
infinity. This is a satisfactory approximation if the dispersions are small compared to the 
mean values, which is one of the basic requirements mentioned above. 
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1=1 



1 + 



TV 



cxp 



/ 2 q,-J /x CJjq,/(2T 2 ) 



18) 



Thus, for the approximation (Q) to be valid, it is necessary that CJ^/(TV) <C 1 and 
J M CJ^/T 2 <C 1. The Jacobian matrix J M can be calculated by taking the derivative of 
the quantum numbers conservation constraint: 

d 



E ^3 n i = E <i? E n 3i ^ E <?j 



dfi k 



80° ^ 3 3 

u ^Ci j j 1=1 ^ k 



(19) 



where is the i unitary vector. If we define the matrix B: 



j i=i 



(20) 



then the righthand equality in eq. (|I9|) can be inverted so to obtain the derivatives of the 
chemical potentials: 



T(B 



(21) 



It should be noted that the matrix B would be equal to 2A if [i jT = 0. Since fx/T is 
generally 0(1) it turns out that B _1 = (9(A _1 ), hence it is expected to be much smaller 
than 1. To complete the Jacobian matrix J M , we take the derivative of the quantum 
numbers conservation constraint with respect to V for V = V, yielding: 



Bdfi Q° 1 „ - 2 „ 

— — h -=- = > > t q,-A L qjnu = 

Tdv v w y fri Hj ^ 3 



(22) 



We use the Boltzmann limit for all hadrons in the last term, namely we keep only the 
first term of the series. By using this approximation, which is satisfactory if T ~ 170 
MeV, we conclude that the last term is ~ (1/4)B -1 Q°/V which is much less than Q°/V. 
Therefore: 



^ c -TB-2! 
dV V 

Finally, the Jacobian matrix J M turns out to be: 



(23) 



J ~ 

-V — 



T(C(A" 1 ), O(A^), C^A^Q /^)) 



(24) 



where each term is meant to be a column vector. This result can be used in conjunction 
with the eq. fll8[) to establish the validity of the approximation (f|). If A -1 <C 1 moderate 
fluctuations of B, Q, V in comparison with the mean values are needed in order that 
CJ£/TF< 1 and J M CJ£/T 2 < 1. 
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C Fitting Procedure 



We adopted a two-step fit procedure to also take into account the uncertainties on input 
parameters such as hadron masses, widths and branching ratios, which in principle can 
play a significant role in the test of the model. Firstly a x 2 with only experimental errors 
has been minimized and preliminary best-fit model parameters T , Vq, ~p B0 have been 
determined: 

X 2 = U=AVl 2 Vi ] , (25) 
°i 

where the index i runs over the M data points. Keeping the preliminary model parameters 
fixed, the variations Ay' theo of the multiplicities (or ratios) corresponding to the variations 
of the I th input parameter by one standard deviation have been calculated. Such variations 
have been considered as additional systematic uncertainties on the multiplicities and the 
following covariance matrix has been formed: 

r;f Y^^/M (26) 

I 

to be added to the experimental covariance matrix C exp . Finally a new x 2 has been 
minimized: 

At 

x 2 = E (yT p - <r°)[(c exp + en -.ud/r - yf co ) ( 27 ) 

from which the best-fit estimates of the model parameters and their errors have been 
obtained. Actually more than 130 among the most significant or worst known input 
parameters have been considered and the corresponding la variations performed. This fit 
technique upgrades the one used by one of the authors in the analysis of thermal hadron 
production in e + +e~, p + p and p + p collisions [|, [7|, H in that the off-diagonal elements 
of C sys are also included. 



D x 2 an d Data Redundancy 

We prove that the parameters fitted with a x 2 minimization and their errors are not 
affected by the presence of redundant data. Let y\ . . . be a set of experimental mea- 
surements among which yk- ■ - Vn are measurements of the same variable. Let y = f(x, a) 
be the functional dependence to be tested where a is a set of parameters to be determined 
by means of a x 2 minimization: 

x2 = ^ (m-f(x l ,*)) 2 (2g) 



i=i °t 



The eq. (EH) can be written also: 



2 %^ (Vi~ f(xi,&)) 2 (y l -y + y- f(xk,a)) 2 
X - 2^ Z2 h 72 ' 



i=i °i i=k a i 
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where y is the weighted average of yk ■ ■ ■ Vn', ah these values correspond to the same 
abscissa x^. Hence: 



2 ^ (yi-M.a)) 2 t ^ { yi -yf ^ {y-f{x k ^)Y in ^ { yi -y)(y-f(x k , & )) 
X + ^^^ + -^ + • 

i=l u « i=k U i i=k U i i=k u t 

(30) 

The second term in the above equation is simply the x 2 of the weighted average while the 
third term can be written (y — f(x k , a)) 2 /<r|, Oy being the error on the weighted average 
y; the fourth term vanishes by definition of weighted average. Therefore: 

X 2 = Xwa + X% , (31) 
where Xwa * s the x 2 of the weighted average and: 

2 t, 1 (?/ t -/(x,,a)) 2 (y- /Qr fc ,a)) 2 

X/it = L T2 + 3 ( 32 ) 

i=l U « U J/ 

is just the correct x 2 to minimize, for the N — k + 1 redundant points have been replaced 
with their weighted average. Since Xwa does not depend on a, the minimization of either 
X 2 or Xfiti the latter being the correct one, leads to the same results. On the other hand, 
if n — dim(a) is the number of fitted parameters, the normalized x}u * s: 

2 X ~~ Xwa / Q q\ 

x * = ~T^r (33) 

instead of x 2 /(iV — n). 
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Table 1: Comparison between fitted and measured hadron abundances and ratios. All 
quoted multiplicities do not include feeding from weak decays unless otherwise stated. 
Note: the x 2 's calculated by using values quoted below differ from those of Table ^ as the 
latter include contribution from uncertainties on input hadron parameters. 



Hadron 


Measured 


Fitted 


Fitted 
with 7s =1 


Reference 


S+S collisions 



h" ( a ) 


98±3 


K+ 


12.5±0.4 


K" 


6.9±0.4 


K° 


10.5±1.7 


A( b ) 


9.4±1.0 


A( b ) 


2.2±0.4 


P-P ( c ) 


21.2±1.3 


P( d ) 


1.15±0.4 



92.63 
12.68 
7.611 
9.939 
7.692 
1.474 
21.49 
2.092 



82.04 
13.75 
7.785 
10.49 
10.13 
2.825 
19.79 
2.314 




B 



ff2 



B 
B 
[0 
B 



S+Ag collisions 



A( b ) 
A( b ) 

P-P ( c 
P( d ) 



186±11 
15.5±1.5 
15.2±1.2 
2.6±0.3 

43±3 
2.0±0.8 



171.3 
17.43 
13.99 
2.223 
43.44 
3.381 



147.2 
19.44 
17.44 
2.612 
39.18 
2.401 



B 

12 



B 
B 
B 
B 



Pb+Pb collisions 



Net baryon 


372±10 


375.7 


372.6 


h" ( a ) 


680±50 


650.2 


638.5 


K° 


68±10 


58.27 


73.44 


4> 


5.4±0.7 


5.759 


5.648 


p-p ( c ) 


155±20 


155.3 


147.8 


K+/K" 


1.8±0.1 


1.652 


1.700 


A/A 


0.2±0.04 


0.188 


0.016 



B 
U 



[0 

B 



a - Defined as tt" + K~ + p 
b - Includes feeding from H 

c - Measured with the '+ - — ' method, in this case limited 
rapidity acceptance (0.2-5.8) to exclude spectators 
d - Measured in a restricted rapidity interval and extrapolated 
by assuming that p has the same rapidity distribution as the A 
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Table 2: Hadron gas model fitted parameters. The first set of parameters has been ob- 
tained with a three-parameter fit by setting 73 = 1. The second set is the four-parameter 
fit result when only experimental errors are used (first step of the fitting procedure) while 
the last set is the final result including uncertainties on masses, widths and branching 
ratios. Also quoted are the obtained chemical potentials, the matrix A -1 elements and 
the x 2 - The x 2 f° r S+S collisions within brackets is its corrected estimate accounting for 
kaons data redundancy (see text). 



Parameter 


S+S 


S+Ag 


Pb+Pb 


T (MeV) 
FT 3 exp[-0.7GeV/T] 
7s (fixed) 

~Pb/ t 
X 2 /dof 


208.3±10.4 
2.782±0.091 
1 

1.145±0.066 
34.0/5 


179.9±7.8 
4.91±0.30 
1 

1.470±0.080 
22.3/3 


125.4±4.6 
13.1±1.1 
1 

2.404±0.14 
22.5/4 


T (MeV) 
\/ T 3 exp[-0.7GeV/T ] 

7so 
XoVdof 


182.4±9.2 
3.51±0.15 
0.732±0.038 
1.248±0.074 
17.1/4 


181.8±6.9 
6.20±0.45 
0.727±0.057 
1.365±0.072 
7.74/2 


192.6±8.1 
24.3±1.6 
0.616±0.043 
1.207±0.071 
3.99/3 


T (MeV) 
FT 3 exp[-0.7GeV/T] 

7s 

tVt 

X 2 /dof 

7%/T 
¥q/T 
A- 1 


180.5±10.9 
3.48±0.16 
0.747±0.048 
1.22±0.10 
12.4/4 (11.9/3) 
-0.320 
-0.00217 
(-2.00 + 6.25)10" 2 


178.9±8.1 
6.29 ±0.47 
0.711±0.063 
1.350±0.081 
6.44/2 
-0.363 
-0.0316 
(-1.10 + 3.54)10" 2 


192.6±19.3 
24.3±2.2 
0.620±0.049 
1.21±0.12 
3.16/3 
-0.372 
-0.0655 
(-0.30H-0.87)10- 2 
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Figure 1: Particle ratios in the grand-canonical approximation for 73 = 1. The bands 
correspond to ±ler deviations of the experimental ratios summarized in Table [I]. The 
symbol K stands for (K) = (K+) + (K~) + 2(K°). 
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Figure 2: Intensive thermal fit parameters as a function of system size. Also plotted T 
and 7s fitted in p + p collisions at y/s = 19.5 GeV ||. 
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Figure 3: The strangeness suppression factor As = 2(ss)/((uu) + (dd)) in high energy 
collisions as a function of centre of mass energy (nucleon-nucleon centre of mass energy for 
heavy ion collisions) calculated within the off-equilibrium hadron gas model. For e + + e~, 
p + p and p + p collisions the ratios have been calculated by using model parameters 
quoted in Ref. [p7 |. For p + p in order to estimate the possible influence of the annihilation 



process, we plotted in addition the As value calculated by including initial valence quarks 
and antiquarks (lower points). For e + + e~collisions the leading s quarks in e + e~ — > ss 
have been subtracted to calculate As. 
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Figure 4: The strangeness suppression factor As = 2(ss)/((uu) + (dd)) in a hadron gas at 
fixed T, V and 7s as a function of the number of initial protons (baryon number equal to 
electric charge). 
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